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Abstract 

Despite the intense theoretical and experimental effort, an understanding of the superconduct- 
ing pairing mechanism of the high-temperature superconductors, leading to an unprecedented high 
transition temperature T c , is still lacking. An additional puzzle is the unknown connection between 
the superconducting gap and the so-called pseudogap which is a central property of the most un- 
usual normal state. Starting from the t-J model, we present a microscopic approach to the physical 
properties of the superconducting phase at moderate hole-doping in the framework of a novel renor- 
malization scheme, called PRM. This approach is based on a stepwise elimination of high-energy 
transitions using unitary transformations. We arrive at a renormalized 'free' Hamiltonian for the 
superconducting state. Our microscopic approach allows us to explain the experimental findings 
in the underdoped as well as in the optimal hole doping regime. In good agreement with exper- 
iments, we find no superconducting solutions for very small hole doping. In the superconducting 
phase, the order parameter turns out to have ci-wave symmetry with a coherence length of a few 
lattice constants. The spectral function, obtained from angle-resolved photoemission spectroscopy 
(ARPES) along the Fermi surface, is also in good agreement with experiment: The spectra display 
peak-like structures which are caused alone by coherent excitations in a small range around the 
Fermi energy. 

PACS numbers: 71.10.Fd, 71.30.+h 
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I. INTRODUCTION 



Since the discovery of superconductivity in the cuprates 1 , enormous theoretical and ex- 
perimental effort has been made to investigate the superconducting pairing mechanism which 
leads to an unprecedented high transition temperature T ( £Ji. The generic phase diagram of 
the cuprates shows a wide variety of different behavior as a function of temperature and level 
of hole doping. In particular, with increasing hole doping away from half-filling, the physical 
properties completely change at the transition to the superconducting phase. A large num- 
ber of experiments using angle-resolved photoemission spectroscopy (ARPES) have revealed 
a strong momentum dependence of the superconducting gap^~i2. An additional puzzle is the 
unknown connection between the superconducting gap of the superconducting phase and 
the so-called pseudogap which is a central property of the most unusual normal state of the 
cuprates. 

Superconductivity is usually understood as an instability from a non-superconducting 
state. Therefore, often in theoretical investigations, the starting point was either the Fermi- 
liquid or the anti-ferromagnetic phase at large or low doping. In this paper, we take a 
different approach and only consider hole fillings, in which either a superconducting or a 
pseudogap phase is present. A generally accepted model for the cuprates is the t-J model 
which describes the electronic degrees of freedom in the copper-oxide planes for low ener- 
gies. Alternatively, one could also start from a one-band Hubbard Hamiltonian as a minimal 
model. However, for low energy excitations, the latter model reduces to the t-J model, so 
that both models are equivalent. In a preceding paper—, henceforth denoted by I, we have 
investigated the pseudogap phase in the cuprates on the basis of the t-J model. Our aim is to 
extent the microscopic approach from paper I to the superconducting phase. As our theoret- 
ical tool, we use a recently developed projector-based renormalization method which is called 
PRM^. The approach is based on a stepwise elimination of high-energy transitions using 
unitary transformations. We thus arrive at a renormalized 'free' Hamiltonian for correlated 
electrons which can describe both the superconducting phase and the pseudogap phase. For 
the superconducting phase, the order parameter turns out to have <i-wave symmetry with 
a coherence length of a few lattice constants. The basic feature for the understanding of 
the superconducting pairing mechanism in the underdoped regime is a characteristic elec- 
tronic oscillation behavior between neighboring lattice sites. The oscillation becomes less 
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important for larger 5 which agrees with the weakening of the superconducting phase for 
larger hole doping. The spectral function, obtained from angle-resolved photoemission spec- 
troscopy (ARPES) along the Fermi surface, also agrees well with experiment: The spectra 
display peak-like structures which are caused alone by coherent excitations in a small range 
around the Fermi energy. 

After a short introduction of the model in Sec. II, we apply the projector-based renor- 
malization method (PRM) in Sec. Ill to the t-J model. The results will be discussed in 
Sec. IV. 

II. MODEL 

In the preceding paper I, we have investigated the pseudogap phase in the cuprates on 
the basis of the t-J model. We adopt the same model also for the superconducting phase of 
the hole-doped cuprates. As before, we restrict ourselves to moderate hole concentrations 
away from half-filling outside the antiferromagnetic phase 

ij,a ia ij 

The t-J Hamiltonian consists of a conditional hopping term and an antiferromagnetic ex- 
change interaction and acts in a unitary space with empty and singly occupied sites. The 
Hubbard creation and annihilation operators c£? = c^(l — ify &) in Eq. ([T]) obey nontrivial 
anti-commutator relations 

[4, c jV ]+ = <% (<w£U*") + s^s?) . (2) 

Here, S q is the local spin operator and T> a (i) is defined by T> a (i) = 1 — n it _ a . In Fourier 
notation, the t-J model ([I]) reads 

+ £ J q S q S_ q . (3) 

k,a k q 

£k measures the one-particle energy from the Fermi energy = — tije l ^ R ^ Rj ^ — /i. 

Note that in Eq. ([3]), we have introduced an infinitesimal field Ak,A —* which breaks the 
gauge symmetry in the superconducting phase. 
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III. RENORMALIZATION APPROACH FOR THE SUPERCONDUCTING 
PHASE 



Let us apply the PRM to the t-J model in the superconducting phase. We consider 
the case of moderate hole-doping, where superconductivity occurs. As before, the hopping 
element t between nearest neighbors is assumed to be large compared to the exchange 
coupling J. Therefore, we can decompose the Hamiltonian into an 'unperturbed' part H,q 
and into a 'perturbation' Tix, 

T~Lo = ^2ekcl a c ka + ( A k,AC k T a k| + A k,A £ k|C k T) + Yl J q^o(q), 

kcr k q 

Hx = ]TJ q (A(q) + 4(q))- (4) 

q 

The decomposition (jl]) is an extension of the former decomposition for the pseudogap phase 
to the superconducting phase. It is based on a splitting of the exchange into two parts. 
The first one, containing A$, commutes with 7i t and should, therefore, be a part of the 
unperturbed Hamiltonian 7{q. In contrast, the two operators A\ and A\ do not commute 
with Ti t and belong to TL\. They are defined by 

Mi) = \ (s q s-, + -^s q s_ q ) , (5) 
A(q) = i( Sq _i_s q ) (s-,-^s_ q ), 
-4!(q) = j(s q + ^s q ) (s- q + £s_ q ). 

and obey approximately the following relations: 

LoA(q) = 0, L ^i(q) = 2£D q ^i(q), L 4(q) = -2w q 4(q). (6) 

Here, L is the Liouville operator corresponding to TCq, where L is defined by L C = [7i ?C] 
for any operator variable C, and u) q is given by 

<% = 2P (^ - t\) = u\ > 0, t\ = 4 ■ (7) 

J(#) 

A. Renormalization equations 

The derivation of the renormalization equations for the parameters of the Hamiltonian 
runs parallel to that for the pseudogap phase. The aim of the projector-based renormaliza- 
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tion method (PRM) is to eliminate all transitions due to 7i\ between the eigenstates of Tio 
with non-zero transition energies. Let us assume that all excitations with energies larger 
than a given cutoff A have already been eliminated. Then, an ansatz for the renormalized 
Hamiltonian 7i\ should have the following form, 

T^A = ^0,A + 7^i,a (8) 

with 

H ,x = H t ,x + J q> A -V(q) - ( Ak . A £ M 6 -k,i + A k,A c-k,iCk, T ) + E x , (9) 

q k 

Hl >X = ^J q , A 0(A-|2^ q , A |)(A,A(q)+^ t 1 , A (q)) • 
q 

Ti-t,\ — J2\ca £ k,A c^o-CkCT is the renormalized hopping term and depends on A. The other 
parameters A k A , £> qiA , and J qiA in Eq. (JH]) are also A-dependent. However, the A-dependence 
of J qjA can be suppressed according to paper I. 

The A-dependent operators A a> \(ci) (a = 0, ±1) in Eqs. ([9]) are defined as in Eqs. ([5]). 
However, S q and cD q have to be replaced by S q y \ and cD qiA , 

S q ,A = ^[Wo,A,S q , A ] wi^wj, (10) 

a2 _ 9P n f7 2 -t 2 ) t 2 - \^ t 2 P *»( R *- H *) 

^ q ,A — Z -M) l r q =0,A r q ,Aj ■> r q ,A — 2.^1 *i> A 

In order to derive renormalization equations for the parameters of H,\, we eliminate 
all excitations within an additional energy shell between A and a reduced cutoff A — AA. 
According to paper I, this is done by applying a unitary transformation to TC\, 

W(a-aa) = e x ^H x e- x ^. (11) 

The generator X A) aa was constructed in paper I and is given in lowest order perturbation 
theory by Eq. (1.37), 

*A,AA = £ AA ) (AA(q) " 4,A(q)) • (12) 

Here, q (A, AA) denotes a product of two O-functions 

6 q (A,AA) = e(A-|2^ q , A |)0(|2cu q , A „ AA |-(A-AA)), 
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which confines the elimination range to excitations with |2c<j qiA _AA| larger than A — A A and 
[ 2a>q 5 A | smaller than A. Roughly speaking, for the case of a weak A-dependence of |^ q ,A|, the 
elimination is restricted to all transitions within the energy shell between A — AA and A. 
According to Eqs. the generator X x ,a\ can also be expressed by 

X XA X = i^- @ ^ A ' AA ) ( S * S-q,A + S q ,A S_ q ) . (13) 

q 1> X 

The explicit evaluation of the unitary transformation fill I) follows that of paper I. In 
perturbation theory to second order in J q , one finds 

Ha-aa = e x ^ H x e' x ^ = H x °l AX + w£ AA + H x \ x + • ■ ■ , (14) 

where 

^i-AA = Ht,X ~ E ( A M 4/-k,| + K,X C-k,|Ck, T ) + E X , 
k 

Wa-aa = E^/ q A,A(q) + [x XAX ,n t , x ] + ^ j q e(A - |2£ q , A |) (Aa(q) + ^U(q)) , 

q q 

^A-AA = 2 t^ A ' AA ' [^A,AA, ^,a] ] + J<\ \X\AXi A,A(q)] 

q 

+ J q 0(A - |2w q , A |) [ X A , AA , A,A(q) + A{ x (q) ]. (15) 
q 

All expressions agree with those of paper I, except that in 7~t\L AX the new symmetry breaking 
terms appear. The commutators can be evaluated as in paper I. Let us at first investigate the 

(2) 

effect of the second order term H A _ AA . The obtained operator expressions have to be reduced 
in a further factorization approximation to operator terms appearing in 7i x . Thereby, also 
a reduction to operators cLcL k , and c_kj,Ckt has to be included. The final result has to be 
compared with the formal expression for 7^ A -aa ; which corresponds to the expression (JHJ) for 
7i x , when A is replaced by A — AA. According to Appendix A, the following second order 
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renormalizations to and to the order parameter Ak,A are found 




k' 

where we have defined 

nZ L) = <4,W " i £<&A*> (18) 

k' 

as non-local part of the one-particle occupation number per spin direction. An equivalent 
equation also exists for E\-ax- The quantity (S q) AS_ qi A) is a correlation function of the 
time derivatives of S q and was evaluated in paper I. Note that an additional contribution to 
£k,A-AA, proportional to the correlation function (S q - S_ q ), has already been neglected. The 
remaining expectation values in ffl6l) . ffTTj) have to be calculated separately. In principle, 
they should be defined with the A-dependent Hamiltonian 7i\, because the factorization 
approximation was employed for the renormalization step from TC\ to Tix-AX- However, 7i\ 
still contains interactions which prevent a straight evaluation of A-dependent expectation 
values. The best way to circumvent this difficulty is to calculate the expectation values with 
the full Hamiltonian TC instead of with 7i\. In this case, the renormalization equations can 
be solved self-consistently, as it was done in paper I. 

Up to now, the renormalization contributions were evaluated from the second order term 
T^x-ax °f T^a-aa- Inserting £k,A— aa and A^a-aa into Eq. ( fl4l) . we obtain 

Hx-AX = Ht,X-AX - ^2 ( A k,A-AA 4,T^-k,J. + ^*k,X-AX C-k,|Ck,t) + ^A-AA + ^A-AA-(19) 

k 

The first order term "H^-ax nas s ^ m ^° ^ e evaluated. This can be done along the procedure of 
paper I. The final result for the renormalized Hamiltonian 7Ya-aa reads Ha-aa = ^o,a-aa + 
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H^x-ax, with 

Kq,X-AX = Kt,X-AX - ^ ( A k,A-AA Ck,T 5 -k,| + A fc,A-AA C-k,iQk,t J + ^A-AA (20) 

k 

,A— Aa(q) j 

q 

Hi,x-ax = X] ^ q @(A - AA - |c£; q> A-AA|) (^i,x-ax((1) + A\ x _ AX (q) 



q 

The renormalized Hamiltonian 7Ya-aa has the same operator structure as 7i\. Therefore, 
we can formulate a renormalization procedure as follows: We start from the original t-J 
model in the presence of a small gauge symmetry breaking field. The energy cutoff of the 
original model is denoted by A = A. Starting from a guess for the unknown expectation 
values, which enter the renormalization equations ffTB]) and ffTTl) . we proceed by eliminating 
all excitations in steps AA from A = A down to A = 0. Thereby, the parameters of the 
Hamiltonian change in steps according to the renormalization equations ( fl6l) and ( fTTI) . In 
this way, we obtain a final model at A = 0, in which the perturbation 7ii t \ is completely 
integrated out. It reads 



n 



A=0 



£ k,A=0 4r ^ - Yl ( A k,A=0 4,T £ -k4 + A t,X=0 C-k.^k.T 

kcr k 



+ J]j q A,A=o(q) + £A=o. (21) 
q 

Unfortunately, due to the presence of the .Ao-term, the result ff2T]) does not yet allow us to 
recalculate the expectation values, since the eigenvalue problem of H\=o can not be solved. 
Therefore, a further approximation is necessary. It consists of a factorization of the second 
term in 



^q-4().A=o(q) — ( SqS-q + — S q .A=(}S- q> A=0 

q 2 V ^q,A=0 



(22) 



According to Appendix A, we end up with a modified Hamiltonian which will be denoted 
by W«, 

Jn 

kcr k q 



= E ^ ) 4. a k cr - £ (£« 4 , T al kii + A«* e__ ka c k , T ) + E y s q s_ q + 

kcr k q 

(23) 

Here, not only the electron energy £k,A=o but also the order parameter A k ,A=o is modified 



according to 



~(1) 1 3J q . , 2 (NL) 

B k - £k,A=0 - 4^2 l £ k,A=0 - ^k+q,A=0j %+q, CT > 

q q>-^=o 



- T7 ^2 Art q ( £k ' A=0 ~~ £ k+q,A=o) 2 (C-(k+q)| C k+qT ) , (24) 



A (!) _ A 

Ak " Ak ' A -" .v^ 



where ri k ^ is defined in Eq. ( TT8i) . Note that the operator structure of Tt^ 1 ) agrees with 
that of the original t-J model of Eq. ([3]) in the presence of the symmetry breaking field. 
However, the parameters have changed. Most important, the strength of the exchange 
coupling in Eq. fl23|) is decreased by a factor of 1/2. This property allows us to start the 
whole renormalization procedure again. We consider the modified t-J model (|23l) as our new 
initial Hamiltonian (at A = A) which again has to be renormalized. The initial values of the 
new Hamiltonian Ti^ at cutoff A = A are e £\ A^, and JJ2. After the new renormalization 
cycle, the exchange coupling of the renormalized Hamiltonian Ti,^ is again decreased by a 
factor of 1/2, until, after a sufficiently large number of renormalization cycles (n — > oo), the 
exchange completely disappears. Thus, we finally arrive at a 'free' model 

W = Yl §k 4* Ck. - ( Ak + A k C-k,|Ck, T ) + E . (25) 

kcr k 

Here, we have introduced the new notation, Ji = 7^( n ^°°) ) £ k = e^ - * 00 ^ Ak = A k n ^°°- ) , 
and = ^( n "^°°). Note that the Hamiltonian Iri allows us to recalculate the unknown 
expectation values. With these values, the whole renormalization procedure can be started 
again, until, after a sufficiently large number of such overall cycles, the expectation values 
converge. Then, the renormalization equations have been solved self-consistently. However, 
the fully renormalized Hamiltonian f l25l) is actually not a 'free' model. Instead, it is still 
subject to strong electronic correlations which are built in by the presence of the Hubbard 
operators. 



B. Evaluation of expectation values 

The expectation values in Eqs. (Tl6|) . (fTTl) . and (1241) are formed with the full Hamiltonian. 
To evaluate an expectation value (A), we have to apply the unitary transformation also on 
the operator variable A, 

(A) = ^ e _ m } = (A(X)) nx = (A)n , (26) 
9 



where we have defined .4(A) = e Xx Ae Xx and A = A(X — > 0). Thus, additional renormal- 
ization equations for A(X) have to be derived. 

1. ARPES spectral functions 

First, let us consider the spectral function from angle resolved photoemission (ARPES). 
It is defined by 



A(k,u) = 




and can be rewritten by use of the dissipation-fluctuation theorem as 

A(k,u) = ^-L^-QfG^w). (28) 

Here, ^(^(k, uj) is the dissipative part of the anti-commutator Green function, 

1 f°° 

9fG(k,w) = — / ([cl(-t), c ka } + )e^dt = ([cl,5(l + cu)c ka } + ). 

The time dependence and the expectation value are formed with the full Hamiltonian Tl, and 
L is the Liouville operator corresponding to Tl. According to Eq. ( 1261) . the anti-commutator 
Green function can be expressed by 

9fG?(k,w) = <[c| C(T (A),5(L A + ^)c k(J (A)] + > A , (29) 

where the creation and annihilation operators are subject to the unitary transformation. In 
order to derive renormalization equations for Ck CT (A) and c kcr (A), we restrict ourselves to a 
weak coupling theory. In this case, all contributions to the unitary transformation from the 
symmetry breaking fields can be neglected. Therefore, we can take over the previous ansatz 
(1.59) for Cko-(A) from paper I: 

Cka(A) = M k ,ACk CT + 7^7 ^ ^k,q,A ^Tf" ^(<?a/3 ' <?<T 7 )(£k',A ~ ^k'+q.A) C k '+qa ^k'/3 C k +q 7 - 

qk' 1' A a/?7 

(30) 

Note that the dominant A-dependence of Cko-(A) is transfered to the parameters Wk,A an d 
^k.q.A- The general renormalization scheme was already established in paper I. Thus, running 
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through the renormalization cycle many times (n — > oo), the exchange interaction will 
completely be eliminated. For n — > oo, we arrive at the fully renormalized operator 

d£+°°)(X = 0) = W k C k(T + ^ 5^Uk,q ^ • ^<r 7 )(£k' - £k'+q) 4'+qa Ck'/3 Ck+q 7 , 

qk' q a/3-y 

(31) 

where w k = u^^^\ £ k)q = v^Z'£j Q , and = £^^*^ ■ Using the renormalized Hamiltonian 
7i of Eq. (1251) . the spectral function 3G(k,w) can be transformed to 

9G(k, W ) = <[4r OO)t (A = 0),5(L + ^)c(r° o) (A = 0)] + >^ (32) 

where the Liouville operator L is related to 7i. The expectation value has to be eval- 
uated with 7i. For this purpose, we introduce new approximate quasiparticle operators 
(Appendix B), 

a l = U k c k T - V k c_ k4 , 

4 = U k cl k4 + V k c kjT , (33) 

which fulfill the following relations: La k = E^a^ and L/3 k = £/ k /3 k , where i? k = 
-D 2 A k . Inserting Eq. fl3Tl) into Eq. (1321 and replacing all c k( ^-operators by the quasi- 
particle operators a^ 1 and /3 k , the ^-functions can be evaluated. For the expectation values, 
we restrict ourselves to the leading order in the superconducting order parameter. The re- 
sulting expression for SG(k, uj) reads: 

9fG(k,w) = ^ \ (l + 5(u - E k ) + (l - 6(u + E k )\ (34) 



- I V E k J \ /. k 

T ~ \ 2 



3D x ^ 

qq' 



4cj 

q 



o— I (^k+q' — ^k+q+q') 2 



x {n k+q+q /m k+q / + n k+q (£> + n k+q / - n k+q+q /)} 

oT^TT ^k,qSk,q' (^k+q' ~ ^k+q+q') ( £ k+q ~ ^k+q+q' 
X {(™k+q' - ^k+q)^k+q+q' ~ ™k+q' (™k+q + D)} 

x5{u + sign(e r k+q+q /) J B k+q+q . - sign(£ k+q ,)£ k+q / - sign(e k+q ).E' k+q } 
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where and are defined by = (cj^c^)^ and = (ck CT c ko .).£. For and m^, we 
use the Gutzwiller approximation^, 

n k = (D-q)+qf(e k ), (35) 

1 — ?7 

m k = g(l-/(e k )) with g=- — , 

1 — n/2 

where /(e k ) is the Fermi function, /(e k ) = 0(— e k ) for T = 0. Note that mjj is proportional 
to the hole filling 5 — 1 — n. Obviously, the application of c k(T on a Hilbert space vector is 
non-zero only when holes are present. In contrast, does not vanish even at half-filling. 



2. Pair correlation function 

In order to evaluate the superconducting order parameter Ak, we have to know the 
superconducting pairing function (c_ k | c k |). Here, the expectation value is defined with 
the full Hamiltonian for the superconducting phase. We first have to transform the pairing 
function, according to Eq. (1261) 



(c-kiCkt) = (c- k |(A)c kT (A)) 



where the expectation value is now formed with the Hamiltonian Tlx, given by Eq. (jSJ). In 
a weak coupling theory, all contributions from the symmetry breaking fields to the unitary 
transformation of c_kj.(A) and Ckf(A) can again be neglected. Therefore, we can immediately 
take over our previous result (J30l) for c k<T (A). For the full renormalization (n — > oo), we 
obtain 

(c_k|CkT) = Mk(c-k|Ckt)w ( 36 ) 
3 / J \ 2 

+ 2iV2 I j ~ e k'+q) 2 mk'+q^k'(C-(k+ q )|C (k+q ) T )^ . 

qk' ^ 1' 

Contributions from third order in the superconducting order parameter have been neglected. 
The expectation values on the right hand side are formed with the fully renormalized 
Hamiltonian 7i (Eq. fl25|) ). Using again the approximate Bogoliubov transformation of 
Appendix [Bj we find 

<^fi kT >* = w^-iT^y (37) 
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FIG. 1: The superconducting gap function Ak versus k, as obtained from Eq. (|44j) for a square 
lattice with N = 40 x 40 sites. The parameters are 5 = 0.08, t! = 0.4£, T = 0. Note that the gap 
function shows <i-wave symmetry. 

IV. NUMERICAL EVALUATION FOR THE SUPERCONDUCTING STATE 

Superconducting solutions have been obtained by evaluating self-consistently the full 
PRM renormalization scheme for a sufficiently large number of renormalization cycles. We 
have taken the same parameters as for the normal state in subsection V B of paper I, 
t' = OAt, J = 0.2t. 

A. Order parameter 

1. Zero temperature results 

In Fig. CD, the superconducting gap function Ak is plotted in k-space for optimal doping, 
5 = 0.08. In agreement with experiment, the solution shows <i-wave symmetry with nodal 
lines directed along the diagonals of the Brillouin zone from (— 7T, — 7r) to (tt, 7r) and from 
(7T, — 7r) to (— 7r,7r). No s-wave like solutions were found. 
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k x /7t k x /7i 

FIG. 2: The superconducting gap Ak (left panel) and the superconducting pairing function 
(c_kj, Ckf) for the same parameters as in Fig. [1] plotted as a 2d map. 

In Fig. |2l both the superconduction gap function Ak (left panel) and the pair correlation 
function (c_kj, Okf) (right panel) are shown as a 2d-plot for the same parameter values as in 
Fig. [TJ Again, in both functions, the nodal lines are clearly seen. Moreover, the absolute 
value of the pair correlation |(c_kj. Ck|)| has a pronounced maximum along the Fermi surface 
(FS). This behavior can easily be understood from Eq. ( 1371) . For k- values close to the 
FS, k « kp, where < 0(A k ), the quantity |(c_kj.Ckj)| is of order 0{1). In contrast, 
for k-vectors away from the FS (with ^> (9(Ak)), the pair correlation function is of 
order O(Aft). Note that the gap function |Ak| has only a weak minimum at the Fermi 
surface. Additional weak maxima can be detected for the following k-vectors: (±7r, ±0.557r), 
(±0.55vr,±7r), (±0.5vr,0) and (0,±0.5tt). 

Fig. [3] shows the superconducting gap function Ak on the Fermi surface as a function of 
the Fermi surface angle for three doping values, 5 = 0.05 (underdoped case, blue line), 
5 = 0.08 (optimally doped, black line), and 5 = 0.12 (overdoped, red line). The angle was 
already defined in paper I in the inset of Fig. 3. In all three cases, Ak shows a characteristic 
overall increase from the nodal (<j) = 0) to the anti- nodal point. Note, however, that the 
maximum value is already reached at a finite angle of about 27°, which is followed by a weak 
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FIG. 3: Superconducting gap function Ak (in units of 2t) as a function of the Fermi surface angle 
(j) which was defined in the inset of Fig. 3 of paper I for three doping values, 5 = 0.05 (underdoped 
case, blue line), 5 = 0.08 (optimal doping, black line), and 5 = 0.12 (overdoped case, red line). 

decrease of A k . 

According to Fig. [TJ the gap function shows a pronounced k-dependence in the whole 
Brillouin zone. By Fourier transforming Ak to the local space, 

^ = ^E A ^ )eik(Rl ~ Rj) ' ( 38 ) 

one finds the spatial dependence shown in Fig. HI The figure again reveals the <i-wave 
character of the superconducting order parameter. Note that the strong k-dependence of 
Ak maps on a short range behavior in local space. As is clearly seen, the local order 
parameter decays in space within a few lattice constants. This feature is consistent with 
the experimentally found superconducting coherence length in the cuprates of the order of 
a few lattice constants. The order parameter changes its sign by proceeding along the x- or 
2/-axis. This can be seen for various hole fillings in Fig. O where Ay is shown as a function 
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FIG. 4: Superconducting order parameter in local space Ajj (in units of 10~ 2 (2i)) for optimal 
doping 5 = 0.08 and T = 0. The hopping parameter t! between next-nearest neighbors is given by 
t' = OAt. Rfj and R^- denote the x and y components of Rj — Rj. 

of Rfj (for fixed Rfj = 0). Here Rfj and R\j are the components of the difference vector 
Rjj = Ri — Rj between lattice sites Rj and R^. The alternating sign of Ay seems to be 
reminiscent of the sign behavior of antiferromagnetic correlations. However, the sign change 
is a property of the superconducting state and not of antiferromagnetic correlations. 



2. Finite temperature results 

In Fig. [6], the local order parameter A^ is plotted as a function of T for different values 
of the distance between local sites, k = |Ry|. The curves are obtained from Fourier back 
transforming Eq. fl38|) together with the temperature dependent expression for (c_k|C k j) 
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FIG. 5: Superconducting order parameter Ajj (in units of 2t) in local space along the x direction 
(for fl* = 0) for three different hole fillings 5 = 0.05 (blue), 0.08 (black), and 0.12 (red). The 
parameters t' and T are the same as in Fig. [TJ 

from Sec. IIIIBI All curves vanish at the same temperature T/2t ~ 0.026, which defines 
the critical temperature T c . Note that the temperature dependence of Ajj and thus of 
the gap function Ak resembles that of the order parameter in BCS superconductors. This 
property can be traced back to the diagonalization approach on the basis of a Bogoliubov 
transformation in Appendix [B], which is applied to the renormalized Hamiltonian 7i in the 
superconducting state. Also the pair correlation function (c_kj,Ckf) is evaluated in this way 
which results in a temperature dependence as in BCS superconductors as well. 

In Fig. [3, the critical temperature T c is given as a function of the hole doping 5. The 
parameter values are again if = OAt and J = 0.2t. Note that for small hole doping 5 < 0.03, 
no superconducting solutions are found. Also this result of the PRM is in good agreement 
with experiments. In the underdoped region for 5 > 0.03, the critical temperature T c 
first increases substantially until it arrives a maximum value at about 5 ~ 0.08. Above 
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FIG. 6: Local order parameter Ajj(T) as function of T (both in units of 2t) for different values of 
the distance k = |Rjj|. Note that all curves vanish at the same critical temperature T c . 

the optimal doping concentration of 5 = 0.08, the critical temperature decreases again 
(overdoped region). Within the parameter range, given in the figure, the T c behavior agrees 
very well with experiment. For still larger values of 5 (5 > 0.15), our PRM result for T c 
remains finite. This feature is in disagreement with experiments, where the superconducting 
phase vanishes above a critical hole concentration. However, this defect of the present 
approach is by no means surprising. As was discussed in Sec. IHIl we have argued from 
the beginning that the present approach is not applicable for the case of large hole doping. 
Nevertheless, Fig. [7] demonstrates that we are able to explain the experimental findings at 
least in the underdoped and in the optimal doping regime. For the present parameter values, 
the maximum of T c at optimal doping is approximately given by T c pa 0.061 Assuming a 
bare bandwidth of 8t ~ lO^K, this T c -value corresponds to a critical temperature of order 
50 — 100 K, which is in the correct order of magnitude. 
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FIG. 7: Critical temperature T c as a function of the hole doping 5 for t' = OAt and J = 0.2t. 
No superconducting solution is found for 5 < 0.03. This result explains the vanishing of the 
superconducting phase in the cuprates at very low doping. 

3. Discussion 

Next, we want to discuss the origin of the superconducting pairing mechanism. Let us 
start with the superconducting order parameter after the first renormalization step. 
According to Eq. (I24|) . we have 

A^ _1) = A k;A=0 - — ^ q (£ k ,A=0 - £k+q,A=o) 2 (C-(k+q)| C k+qT ). (39) 

q 1)^=0 

The first term on the right hand side results from second order renormalization contributions 
according to Eq. (IT7|) . The numerical evaluation of Eq. fl39l) shows that it is small compared 
to the second term. According to Sec. IHIl the latter one results from the factorization of the 
contribution ~ S q S_ q in the renormalized Hamiltonian H\=o = ^ q (<^q)/(2iiq)S q S_ q + • • • 
after the first renormalization cycle. Therefore, we can conclude from (IA3I) that the dominant 
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part of the microscopic pairing interaction is given by 

n (sc) = ^ E 7^( £k ~ £ k- q ) 2 (cj cT £t -k|C-(k- q )|C k _ qT + 2cj cT cL ki c k _ qi c_ (k _ q)T ) . (40) 

qk 1 

Here, spin-singlet pairing was assumed. The expression ( HU1) is our central result for the 
superconducting pairing mechanism in the cuprates. In contrast to usual BCS supercon- 
ductors, where the pairing interaction between Cooper electrons is mediated by phonons, 
the present result can not be interpreted as an effective interaction of second order in some 
electron-bath coupling. Note that Eq. (j4"U|) results from the part of the exchange which 
commutes with 7i t . An important feature of the pairing interaction is the oscillation fre- 
quency u>q in the denominator of Eq. ( 1401) , 

uj\ = -2P (t q=0 - t\) = £ 2 q > 0, t\ = Y, til cosq(R, - Ri) , (41) 

which enhances the pairing mechanism for small hole doping, since Pq ~ 5. Therefore, the 
pairing interaction is mediated by oscillating hopping processes between nearest neighbors. 
This was discussed in detail in Sec. IV A of paper I. First, an electron hops to a neighboring 
site which is empty. In the second step, it hops back to the first site, since this site was 
certainly empty after the first hop. Thereby, the presence of short range antiferromagnetic 
correlations in the unperturbed Hamiltonian 7i is crucial, since it prevents the hopping to 
more distant sites. 

In order to derive an approximate gap equation, let us again start from Eq. (1391) . When 
we restrict ourselves to a weak coupling theory, the A-dependence of and Uq t \ can be 
neglected: 

= ~jf E - £ k + q) 2 (£-(k +qH Qk+qT ) , (42) 

q 1 

where the first term from Eq. (1391) was already omitted. For a purely qualitative discussion 
of the gap parameter, let us abandon all higher order renormalization effects, which would 
be included in the full renormalization scheme of Sec. IIHI Inserting the former expression 
flU for (c_ ki c kT ) into Eq. (02}, we find 

^ = -4E^K^- £ ^) 2 ^ 2 1 - 2fiE ^ ) A k+q , (43) 

N q 2^/4 +q + Z^ +q 
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where is again given by = y £ k + D 2 A^, and /(-Ek) is the Fermi function /(-Ek) = 
1/(1 + e^ E ^). Moreover, by replacing on the left hand side also A k by Ak, we arrive at the 
following approximate gap equation 

A * * 4^ £ Sf(^ - ul +<l )- 2f ^ A k+q . (44) 

Note that the main features of our numerical results for the full renormalization scheme 
can already be detected from this equation. Due to the doping dependence of u^, shown in 
Fig. 9 of paper I, superconductivity sets in at the same small 5-value, at which becomes 
non-zero. With increasing hole doping, increases, which also leads to a strengthening of 
the coherent excitation in 3G(k, uo). Moreover, superconductivity is favored for low doping 
due to the factor cu q ~ 5 in the denominator of Eq. (jHJ). Both features together, i.e. the 
increase of with 5 and cD q ~ 5 lead to a maximum of T c at a finite doping value which 
is seen in Fig. [7J The property u> q ~ 5 also explains the decrease of T c in the overdoped 
region, since renormalization processes become weaker for larger 5. 

The preference of the PRM to find solutions with c?-wave symmetry for the gap parameter 
can also be understood from the gap equation ( I4"4|) . For an explanation, let us start by 
dividing the sum over q in Eq. (jSj) into two parts with |£k+ q | < |Ak+ q | and |£k+ q | > |A k+q |. 
Omitting the second sum, one finds 



q,|Ek+ q |<|Ak+ q | 2 V £ k+ q + A k+ q 



For most values of k, the neglected sum is smaller by a factor of order A/t. An exception 
are k- values close to the Fermi surface k « kp (with |£ k | < 0(A k )), which will be excluded 
in the following discussion. Here, the sum with |£k+ q | > |Ak+ q | would be larger by a factor 
of order t/A. With respect to Eq. (1451) . those terms of the q sum are most important, 
which have energies |£k+ q | not exceeding |Ak+ q |. For k- values on the diagonal, k x = k y , 
of the Brillouin zone, it can be seen that q-values with q y pa q x ± 7r lead to small energies 
^k+ q ~ and thus to the dominant contributions in Eq. (|45p . Here, the dispersion relation 
£k — — 2t(cos k x a + cos k y a) was used. However, the prefactor J q vanishes in this case. This 
explains the nodal line k x = k y and similarly k x = —k y of the gap parameter in Fig. 
However note that the exchange constant J q changes its sign as a function of q. From this 
behavior, one can conclude that <i-wave symmetry for the order parameter is more favorable 
than s-wave symmetry. 
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FIG. 8: Spectral functions 9G(k, oS) in the superconducting phase at optimal doping, 5 = 0.08 for 
two fixed A^-values: (a) k x = ir (anti- nodal region) and (b) k x = 5ir/8 (in between anti- nodal and 
nodal region). By varying k y , the Fermi surface is crossed. The other parameters are t' = 0.4i, 
J = 0.2i, and T = 0. 



B. ARPES Spectral functions 

Finally, let us discuss the ARPES spectral function in the superconducting phase. This 
quantity is obtained from the dissipative part of the anticommutator Green function (I28p . 

In Figs. M - EH our results for the superconducting phase are given which are obtained 
from the numerical evaluation of Eq. fl34l) . First, in Fig. [HI we have chosen as parameters: 
5 = 0.08 (optimal doping), T = 0, t' = OAt, and J = 0.2t. Two cuts with fixed k x and 
varying k y are shown. Thereby the FS is crossed. In panel (a), where k x = tt, the spectra 
belong to k- values in the anti-nodal region, whereas in (b) k x = 5tt/8. Here, a k-region is 
probed in-between the nodal and the anti- nodal point. The spectra in both panels display 
peak-like structures in a small energy range around to — 0. Note that all structures are 
caused alone by the coherent part of SG(k, uj) (first line in Eq. (|3lj) ). which consists of two 
peaks at the positions uj = ±-Ek- For k- vectors, far away from the FS (top and bottom 
plots in Figs. [HJa) and (b)), a dominating peak at u ~ £k is found, which arises from the 
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FIG. 9: Spectral functions 9G(k, u) as in Fig. [8]for a fixed fe^-value, k x = tt/2 . By varying k y the 
Fermi surface is crossed in the nodal region. 



excitations at ±-E k , depending on the sign of £f By approaching the FS, a secondary peak 
arises at u « —e^. An expansion of the prefactors in Eq. (1341) shows that in each case 
the secondary peak has a smaller weight of order (Ak/^k) 2 - Only for k-values on the FS 
(£k = 0), the two coherent peaks have equal weight. They are separated by an energy 
distance, which is given by the gap parameter (2Z)Ak). Note that the gap size is almost 
the same for the two cases of Fig. [HJ A comparison of both panels of Fig. [HJ also shows 
that the secondary peak is more pronounced in the anti-nodal region than in-between the 
ant i- nodal and nodal region. Furthermore, the overall dispersion of of the primary peak 
is weaker in the anti-nodal region than for the case of intermediate k x -va\ues. With respect 
to the incoherent contributions to 3G(k, ou), note that for optimal doping the overall weight 
of the coherent and of the incoherent excitations are approximately the same. However, the 
incoherent part of the spectrum is spread over a much larger frequency range. Therefore, in 
a small u;-range, close to the Fermi level, the coherent excitations are dominant. 

In Fig. [9j the spectral function is plotted in the nodal region for fixed k x = tt/2 and 
different values of k y . Thereby, again the FS is crossed. Note that neither a secondary peak 
nor a superconducting gap is found in the nodal region. Also, the coherent peak moves 
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FIG. 10: Symmetrized spectral functions 9G(k, ui) for k- values on the FS between the nodal (top) 
and anti-nodal point (bottom) for two temperatures (a) in the superconducting phase at T = 0, 
and (b) in the pseudogap phase at T = 0.05i. The critical temperature is T c = 0.03t (underdoped 
case 5 = 0.05). The other parameters are the same as in Fig. [8j 

almost unchanged through the FS, when k y is varied. 

Finally, Fig.[10]shows the results for the symmetrized spectral functions 3G(k, uj) for two 
different temperatures (a) T = (superconducting phase) and (b) T = 0.05t (pseudogap 
phase). The k- values proceed on the FS between the nodal (top) and the anti-nodal (bottom) 
point. The hole concentration is 5 = 0.05 (underdoped regime) which leads to a critical 
temperature T c = 0.03t. In the spectra at temperature T = 0, one recognizes the opening of 
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a superconducting gap for all k- vectors except at the nodal point. The gap size as a function 
of the Fermi surface angle is given by the blue line in Fig. [31 Similar as before, the peak- 
like structure arises from the coherent excitations in a small cu-range around u = 0. For the 
higher temperature, T = 0.05t (pseudogap phase), the system is in the normal state. On a 
substantial part of the Fermi surface, the spectra now show the typical large spectral weight 
around u = 0, indicating a Fermi arc of gapless excitations. The Fermi arc extents over 
a finite k-range. In contrast to the superconducting case, the spectrum is now dominated 
by the incoherent excitations. In the anti-nodal region, they form the pseudogap around 
uj = (see also paper I). Note that the pseudogap in Fig. [TOT b) is about ten times larger 
than the superconducting gap at T = (for the present hole doping 5 = 0.05). Note that 
for both temperatures, the spectra are in good qualitative agreement with recent ARPES 
measurements^,—,—. 

Let us finally make one remark concerning the linewidth of the coherent peaks. As was 
already mentioned in Sec. V of paper I, from the experimental point of view, we would 
expect a temperature dependent broadening of the coherent peaks which is caused by the 
coupling to other degrees of freedom. Such a broadening was not incorporated in the present 
approach. Note, however, that a broadening of the spectra is also produced by the incoherent 
excitations of $SG(\t,u). In order to include a temperature dependent broadening of the 
coherent excitations, we have added by hand a small linewidth in Fig. [TO], which is taken of 
the order of ksT. 

V. CONCLUSIONS 

In this paper, we have given a microscopic approach to the superconducting phase in 
cuprate systems at moderate hole doping. Thereby, a recently developed projector-based 
renormalization method (PRM) was applied to the t-J model. Our result for the supercon- 
ducting order parameter shows ci-wave symmetry with a coherence length of a few lattice 
constants which is in agreement with experiments. In contrast to usual BCS superconduc- 
tors, where the pairing interaction between Cooper electrons is mediated by phonons, the 
superconducting pairing interaction in the cuprates can not be interpreted as an effective in- 
teraction of second order in some electron-bath coupling. Instead, the main contribution to 
the pairing results from the part of the exchange interaction which commutes with the hop- 
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ping Hamiltonian 7i t . The superconducting state naturally arises from a typical oscillation 
behavior of the correlated electrons between neighboring lattice sites due to the presence 
of spin fluctuations. The theoretical results can explain the experimental findings in the 
underdoped as well as in the optimal doping regime. The obtained value of T c at optimal 
doping has the correct order of magnitude. 
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APPENDIX A: FACTORIZATION APPROXIMATION FOR S q S q 

The aim of this appendix is to simplify the operator product S q S_ q which enters the 
expressions ^ for 7^ ,a and TLi,\. As in paper I, we start from the expression 

M-q = JJv £ £(** • **r) £^ iqRl - e***) 5> m (e-^ - e**~) clc^ mS c h . 

(Al) 

The four-fermion operator on the right hand side can be factorized in two different ways: 
One can either reduce it to operators cj^c^ or to operators cj ccr cL k _ cr and C-k- a Ck, CT - The 
first factorization will lead to a renormalization of e^, whereas the second one renormalizes 
the superconducting order parameter A^. In the factorization, we have to pay attention 
to the fact that the averaged spin operator vanishes (Sj) = outside the antiferromagnetic 
regime. Moreover, all local indices in the four-fermion term of ( lAlj) should be different from 
each other. This follows from the former decomposition of the exchange interaction into 
eigenmodes of L t , where we have implicitly assumed that the operators S q and S_ q do not 
overlap in the local space. Otherwise, the decomposition would be much more involved. 
However, it can be shown that these 'interference' terms only make a minor impact on the 
results. 

(i) For the 'normal' factorization, we find 

S q S_ q | (i) = -— ^(e k - £k-q) 2 ((4-qaCk-qa)iVL) (4A*W > ( A2 ) 
ko- 

26 



where we have defined (cl c k(T ) 



■k a ^ka)NL = c^Cka - 0-/N) J2v cLcA'^- The attached subscript 
NL indicates that the local sites of the operators inside the brackets are different from each 
other. In Eq. flA2j) . we have also neglected an additional c-number quantity, which enters in 
the factorization, and the sums over the spin indices in Eq. (1A1[) have already been carried 
out 

(ii) By assuming spin-singlet pairing, we obtain from Eq. (1A1I) . 

S q S_ q | (ii) = ^E(£k-£k-q) 2 (4T^| a -(k-q)|Ck-qT+2^^^ (A3) 
k 

According to Sec. IV A, the expression (1A3j) leads to the main part of the superconducting 
pair interaction. In a factorization approximation, the two contributions in (1A3I) can be 
combined to 

3 



at 



SqS_ q | (ii) = — ^(c7 k - £ k _ q ) 2 |(c_ (k _ q)i C k _ qT ) C^IkJ + h c ] ■ 



(A4) 



Using Eqs. (1A2I) and (1A4I) together with Eq. (1221) . one is finally led to the renormalization 
result (124|) for and A k to first order in J. 

The above factorization can also be used to derive the renormalization contributions 
(TTB|) . (TT71) to £ kj A— aa and A k; A-AA- Using the expressions (1A2I) and (1A3I) . we can first simplify 
the second order renormalization T~C X _ AX of 7^a-aa- In analogy to the results of Appendix B 



in paper I, we arrive at 



n 



(2) 

A-AA 



3D 



J, 



q \2 



6 q (A,AA) 



^ ^(2e k - £ k +q - £k-q)(c k aCka) 
ko 



Sq ' S q 



q 1 



+ (S q ■ S_ q ) — ^(2e k - £ k+q - £ k _ q ) C^C^ J 
ko / 

kcr 

Y X^ 2£k ' ~~ £k ' +q ~~ £ k'-q)(c k 'o-'Ck'o') 
k'o' 

E(^k - £k-q) 2 ((c k -qaCk-qQ)7VL) ( ^j^o- ^ko NL 
ko 



qo 

X 



X 



qo 

X 



X 



- £k-q) 2 ((Ck-q|C k _ qT )) (c^C^ + /l.C.). 



(A5) 
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From (1A5j) . the second order renormalizations to £ k A-AA and A^^-aa can immediately be 
deduced. 



APPENDIX B: BOGOLIUBOV TRANSFORMATION FOR THE SUPERCON- 
DUCTING HAMILTONIAN H 

The aim of this appendix is to diagonalize the renormalized Hamiltonian 7i for the su- 
perconducting phase. According to Eq. (1251) . the Hamiltonian H reads 



n 



^ 4, c kCT - ( Ak a M 6 -k,| + A k C-k,|Ck, T ) + E. (Bl) 

ka k 

Due to the presence of the Hubbard operators in Eq. fIBll) . the usual Bogoliubov trans- 
formation can only be applied approximately. Let us start by introducing new fermion 
operators, 

4 = U k c£ )T -V k c_ ka , (B2) 
Pi = U k cl k4 + V k c kiT . 

We require that a k and are eigenmodes of 7i, 

Lal = E k al, L& = E^. (B3) 

In order to find equations for U k and V k , let us insert the expression (IB2I) for a k into the 



first equation of (1B3j) . 

U k Lc k , T - V k Lc- k4 = £ k (U k 4,-r - V k c_ k)i ) . (B4) 

The two commutators on the left hand side of Eq. (1B4[) will be evaluated separately. For 
the first one, Lc ko - = [ft, c k J , we obtain: 

Lc k T = L<c k ,T ~ E A k'[c-k',|C k ',T,c k ,T]- 

k' 

Here, the Liouville operator L t corresponds to the commutator with the hopping Hamiltonian 
Tt t = ^k c k(T Ckcr, which agrees with the fully renormalized Hamiltonian H in the normal 
state investigated in paper I. Therefore, we can use Uc k * = £ k c kcr and find using the anti- 
commutator relation (j2j) 

L4,T = £ ~4, T - 4w E A L ( e " 4kRj ^T(j) Q,l - e" ikR ^rc JiT ) . (B5) 



28 



AT . 



The quantity A*^ is denned by A* itj = i J] k A*e ik ( Rl - R ^, and V a {j) = 1 -n,--* = V + h ia 
was already given in Eq. (j2j). The main contribution to the second term in Eq. ( ]B5j) is 
caused by the following process: First, two holes are generated at sites i and j before the 
hole at j is annihilated again by a local creation operator in cL.. The arising local projector 
V a {i) will be approximated by its average D = (V a (j)) = 1 — {nj_ a ). Thus, we obtain 

D 



e k cl A - DA*c_ ka , (B6) 

where A* - was Fourier back transformed to A k . A corresponding contribution from the 
last term in Eq. ( 1B5I) vanishes, since (S^) = outside the antiferromagnetic regime. The 
evaluation of the second commutator in Eq. (1B4I) can be done in analogy to the result (1B6I) . 

Lc-k,| = -£kC-u,i - DA k 4 jT . (B7) 

Inserting Eqs. (1B6[) and fjBTj) into Eq. (1B4|) leads to the following two equations for U k and 
V k : 

U k (e k - E k ) + V k DA k = 0, 
-U k DA* + V k (e k + £ k ) = 0. (B8) 

The eigenvalue E k for this system of equations is easily obtained: 



£ k =^~ 2 +D 2 |A k | 2 . (B9) 

The expectation value (c k |C_ k j)^, formed with the superconducting Hamiltonian H, is 
found by solving (1B2j) for and c_ k) j. Using the property (1B3j) . one finds 

4 6 t \ = ^(i — 2 \ 

k ' T - k 'Vw 2£ k V 1 + e^ / V ; 
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